Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 18 May 2010 (MN WI^ style file v2.2) 

Direction Dependence in Supernova Data: Constraining 
Isotropy 



o 



r- 



O 

u 

o 



> 

oo 

00 
(N 

o 
o 



% 



Shashikant Gupta ^'^'^'^ and Tarun Deep Saini ^'^ 

^ Indian Institute of Science, Bangalore, Karnataka, India, 560 012 

^ Raman Research Institute, Sadashivanagar, Bangalore, Karnataka, India, 560 080 

^ shashikant@physics.iisc. ernet. in 

* shashikant&rri.res.in 

^ tarun@physics. Use. ernet. in 



18 May 2010 



ABSTRACT 

We re vise and extend the extreme value statistic, introduced in lGupta. Saini &: Laskaii 
(|2008D . to study directional dependence in the high redshift supernova data; arising 
either from departures from the cosmological principle or due to direction dependent 
statistical systematics in the data. We introduce a likelihood function that analytically 
marginalises over the Hubble constant, and use it to extend our previous statistic. 
We also introduce a new statistic that is sensitive to direction dependence arising 
from living off-centre inside a large void, as well as previously mentioned reasons for 
anisotropy. We show that for large data sets this statistic has a limiting form that 
can be com puted anal ytically. We apply our statistics to the gold data sets from 
iRiess et all (^004) and iRiess et al.l (|2007f ). as in our previous work. Our revision and 



extension of previous statistic shows that 1) the effect of marginalsing over Hubble 
constant instead of using its best fit value has only a marginal effect on our results. 
However, correction of errors in our previous work reduce the level of non-Gaussianity 
in the 2004 gold data that was found in our earlier work. The revised results for the 
2007 gold data show that the data is consistent with isotropy and Gaussianity. Our 
second statistic confirms these results. 
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1 INTRODUCTION 

Observations of high-redshift supernovae of Ty pe la 
l|Perlmutter et al.lll999l : lRiess et al.iri998l . 12001 l20o3 i. along 
with observations of the cosm ic microwave background 
l|Benoit et al.l |2002| : iPage et all [20061 indicate a flat uni- 
verse with an accelerating expansion. The existing evi- 
dence is consistent with a universe that is dominated to- 
day by a cosmological constant term in the Einstein's 
theory of gravity, however, due to fine tuning required 
in this mo del, several alternat i ve explanations have been 
suggested (ISilvestri fc TroddenI l2009l : iFrieman et all l2008l : 
ISahni fc Starobinskvll2006l ). If we phenomenologically treat 
this unknown energy component as an ideal fiuid — called 
the dark energy — then the pressure of this fiuid has to 
be negative, that is, for an equation of state p = wp, we 
have p > and w < 0. In this description the cosmological 
constant model has an equation of state given by p = —p, 
and thus we expect that a phenomenological modeling of the 
cosmological data with a constant w should give lu ~ — 1, 



which turns out to be the case; the current data indicating 
one-third of the total density in the form of dark and bary- 
onic matter and two-thirds in the form of the cosmological 
constant — a model that is known as the ACDM model. 
However, it should be noted that the equation of state for 
the dark energy could differ from ui = — 1 in the past when 
the dark energy is expected to be subdominant (for consis- 
tency with constraints arising from the observations of the 
microwave background). Indeed, in many plausible models 
of dark energy the equation of state approaches w = — 1 
only in the recent past and thus may show only tiny depar- 
tures from TO = —1 at low redshifts, while at high redshifts 
it might remain subdominant and therefore would have a 
weaker observational signature. Although it is perhaps hard 
to extend supernova (SN) observations beyond z ~ 3 or so, 
neutral hydrogen observations of the post-reionization epoch 
may be able to provide constraints on dark energy all the 
way up to 2 ~ 6 IjBharadwai et al.ll2009l ). 

To resolve this issue we require data that is precise 



2 Shashikant Gupta and Tarun Deep Saini 



enough to discern tiny variations in the dark energy. It is 
also required that data be available at a large number of red- 
shifts to constrain the detailed temporal behaviour of dark 
energy. At present the only data that comes reasonably close 
to these requirements is provided by the observations of the 
high redshift type la supernovae (SNe), which are believed 
to be standard candles — a claim that may be doubtful 
considering the relatively poorly understood physics, and 
the possibility of physical mechanisms such as dust in the 
inter-galactic medium that systematically dims them. An- 
other concern is that the SN data sets are usually collated 
from several different experiments that might have slightly 
different systematics, due either to instrumental effects or 
the fact that they observe different directions in the sky; for 
example systematic errors in correcting Galactic or source 
galaxy dust extinction might leave residual anisotropies in 
data. These considerations imply that to obtain precise in- 
formation about the behaviour of dark energy we should first 
have a good knowledge of the statistical properties of SNe, 
both random as well as systematic. 

The standard m odel of cosmol ogy is based on the Cos- 
mological Principle IjPeeblea 119931 ). which states that the 
Universe is homogeneous and isotropic on large scales. Re- 
cent work provides some evidence for what is known as 
the Hubble Bubble IjZehavi et al.lll998l : [jha et al.ll2007f ) in- 
dicating that we might be living inside a large void that 
has a different value of the Hubble constant inside from 
what is outside the bubble. There is eviden ce for such 
large scale voids in the CMB maps as well |Cruz et al.l 



large scale voids m tne UMc maps as wen l^Jruz et ai.l 
20051 : Ide Oliveira-Costa fc Tegmark|]2006l : ICruz et al.ll2006l . 



20071 ). suggesting that such large voids are not implau- 



sible. There is no reason for us to believe that we are 
living at the precise centre of such a void. Therefore, 
if we happen to be living in such a void then vari- 
ous distance measures, such as the luminosity distance 
and the angular diame t er di stance, may not be isotropic. 
iGupta. Saini fc Laskarl (|2008l ) (hereafter GSL08) used the 
extreme v alue statistics to show that the two SN data sets, 
iRiess et a l. (2004) (GD04) and Ries s et al.l l|200'it ) (GD07), 
do show some evidence for direction dependence. Several 
other works have also indicated either systematic problems 
with the high-redshift SN data or directional dependence 
in the SN data and other probes (Kolatt & Lahav 2001'; 



Nesseris fc Perivolaropoulos 2004, 2005a , b; ;Jain fc Ralstoni 
2006al lbl: ICooke fc Lvnden-Belj|20ld ). 



We recently found that our analysis in GSL08 contained 
mistakes due to a coding error. In this paper our main task 
is 1) to provide revised results and to extend the previous 
work by introducing a statistic that marginalises over the 
Hubble parameter instead of using its best fit value 2) to 
introduce a new statistic that has a greater sensitivity to 
the signatures of living off-centre inside a large void. We 
also give asymptotic form for this statistic, which makes it 
easier to use it for large data sets. 

The plan of the paper is as follows. In § 2 we introduce 
the likelihood function marginalised over the Hubble con- 
stant and provide a comparison with what we obtain using 
the Bayesian statistic. In § 3 we introduce the two statistics. 
In § 4 we present our results and § 5 we end with conclusions. 



2 MARGINALIZATION OVER Ho AND M 

For a given SN the measured quantity, the distance modulus 
fj., is the difference between the apparent and the absolute 
magnitude 



fi{z) = m{z) - M , 



(1) 



where the apparent magnitude m{z) is a function of the 
intrinsic luminosity of a SN, the redshift z and the cosmo- 
logical parameters; and M is the absolute magnitude of a 
type la SN. The distance modulus can be expressed in terms 
of the luminosity distance Dl as 



/i(z) = 5 log [Dl{z) /Mpc] + 25 , 



where the luminosity distance is given by 
c(l -I- z) r dx 
Ho Jo h{^) 



Dl{z) 



(2) 



(3) 



where h[z;Q.M,^x) ~ H{z;Q.m,^x) / Ho, and thus depends 
only on the cosmological parameters; the matter density Qm 
and the dark energy density ilx. The prescription for the 
variation of dark energy with redshift has to be specified 
separately. For example, in the ACDM model the energy 
density in the dark energy Oa remains a constant. In Eq[T] 
the dependence of the measured quantity /x on M is linear. 
Since ^ depends on the logarithm of the luminosity distance 
it is clear that /i depends linearly on the logarithm of the 
Hubble parameter Ho. Usually the data are given in terms of 
Eq[2j where the constant M has already been marginalised 
over. Thus, instead of two nuisance parameters we are left 
with only one — the Hubble constant Ho- 

The Hubble parameter could in principle be measured 
using the observations of the low-redshift distance-redshift 
relationship, however, the quantities of real interest are the 
cosmological parameter Hm and Jlx. Thus it would be use- 
ful to marginalise over the nuisance parameter Hq. Although 
this can be done numerically while estimating the cosmolog- 
ical parameters, it turns out that our statistic in § 3.2 does 
not allow this. It is clear that the Hubble constant term can 
be eliminated by considering the difference of two magni- 
tudes 



Mi — Mi = 5 log 



(1 -I- Zi) J^^ dx/h{x) 



{l + Zi)j;'dx/h{x) 



(4) 



If the error on /ii is Gi then the standard error for /ii — /ij 
is given by ct = x'^'i +'^f' where we have assumed that 
the errors on the two magnitudes are statistically indepen- 
dent. If we have a large number, A'', of SNe in the sample 
then dividing it into two equal halves we can form A'^/2 such 
differences (assuming an even A''), and estimate parameters 
that would be independent of Ho- However, the degrees of 
freedom reduce by half in this process, thereby reducing the 
amount of information that can be extracted from data. An- 
other reason for our not choosing this method is the fact 
that we are interested in quantifying direction dependence 
in data and therefore our statistics depend on the direction 
each SN is observed. Specifically, we consider statistics that 
depend on the difference of quantities computed on the op- 
posite hemispheres and then maximizing this difference over 
all the directions. It is clear that such a procedure does not 
allow the use of above method. 

For marginalisation we instead use a method that is 
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based on subtracting the magnitude of an arbitrarily chosen 
low-redshift (so its magnitude depends only on the Hubble 
constant) anchor SN from the magnitudes of SNe in a data 
set and then marginalising over magnitude of anchor SN. 
The resulting likelihood function is derived in the Appendix. 
This is equivalent to marginalising over the nuisance pa- 
rameter Ho with a Gaussian prior centered around its value 
derived from the anchor SN, alongwith the corresponding 
standard deviation. The method can be easily generalized 
to a case where the priors on Hq are specified separately as 
described in the Appendix. 



2.1 Comparison with Bayesian Marginalization 

We now compare the analytically marginalised likelihood 
function in Ea lASI to the results of Bayesian marginalisation 
over Ho- For this purpose we consider parameter estimation 
for GD04. We assume a flat ACDM universe for this exercise. 
The dimensionless Hubble parameter is given by 



h{z) = VnM(l + 2)3+nA: 



(5) 



and flatness implies Qa = 1 — ^m- Thus, the only free param- 
eters are Ho and JIm. The normalized likelihood function is 
given by 



L{x.\D,M,Ho) = 



: exp 



2 



(6) 



(27r)^/2\/detS 

Here Xi — /it ~ fithcory {zi; Ho , i^m) , where the subscript 
i = 1, . . . , N and E^^ — Sij a^ is the covariance matrix. The 
superscript 'T' denotes the matrix operation of taking the 
transpose of a matrix. The posterior probability for param- 
eters SIm and Ho is given by 

L{^\nM,Ho)PiHo,nM) 



P{nM,Ho\fl^) 



P(X) 



(7) 



We choose a uniform prior for Hm in the range ^ (^m ^ 1 
and for Ho = lOO^kms"^ Mpc"^ in the range 0.4 < /i ^ 
1. The marginalised probability distribution for the matter 
density is given by 



P(nM|MO« / L{x\nM,Ho)PiHo,nM)dHo. 



(8) 



Ho 



The probability density is normalized after carrying out the 
integration. A similar probability density function for the 
matter density can be obtained by the Bayesian inversion of 
the marginalised likelihood function given in Eq. A5. Fig- 
ure [T] plots a comparison between the two probability den- 
sities and shows that the two distributions are nearly iden- 
tical. 



3 TWO STATISTICS 

For completeness some details are repeated here from 
GSL08. For our analysis we have considered a flat ACDM 
universe, which can be easily generalized to a more general 
model of dark energy. 
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Figure 1. Comparison of numerically marginalised probability 
density for the Hm with the analytic marginalisation according 
to the likelihood function in Eg lASI The two are indistinguishable. 
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Theoretical 



0.5 



1.5 



^2 

Figure 2. A comparison of theoretical and bootstrap probability 
distributions for simulated data. The data comprises 157 SNe, 
whose positions on the sky were generated randomly. There was 
a mistake in a similar figure in GSL08 where the discrepancy 
between these distributions was larger. 



3.1 



A-,^2 statistic 



We consider subsets of the full data set to construct 
our statistic comprising A^subsct data points. For analytic 
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marginalisation we keep two anchor SNe aside. Since the 
ACDM model fits the gold data sets GD04 and GD07 
well, we first obtain the best fit to the full gold data 
sets using the marginalised likelihood function given in 
Eq IA5I and then for each SN we calculate the residuals 

X^ = [m. - ^''^^^(^M «m)] - [mo - ^''^^'"(^O; Om)], which 

is free of Hq by virtue of Eq|4l The standard error for a SN 
at redshift Zi is 0-^(2^), and we assume that the errors on 
SNe are statistically uncorrelated. 

We define Xr = XMZ-'Vsubsct, where the marginalised 
Xm is defined in Eq. IA7I Xr indicates the statistical scat- 
ter of the subset from the best fit ACDM model. Its ex- 
pectation value is unity (see Appendix for a proof), that 
is (Xr) ~ i- We divide the data into two hemispheres la- 
beled by the direction vector h, and take the difference of 
the Xfl computed for the two hemispheres separately to ob- 
tain Axfi = Xri ~ Xji2i where label '1' corresponds to that 
hemisphere towards which the direction vector h points and 
label '2' refers to the opposite hemisphere. We take the ab- 
solute value of Axrf since we are interested in the largest 
magnitude of this quantity. We then vary the direction h 
across the sky to obtain the maximum absolute difference 

A^2 =max{|Axft|} . (9) 

As shown in GSL08, the distribution of A^2 follows a 
simple, two parameter Gumbel dis tribution, characteristic 
of extreme val ue distribution type I IjKendall fc Stuartll 19771 : 
iGumbell [19651 ). 



Pi\ 



■ exp 



A. 



exp 



' exp 



A^ 



(10) 
where the position parameter m and the scale parameter 
s completely determine the distribution. To quantify depar- 
tures from isotropy we need to know the theoretical distribu- 
tion, which is calculated numerically by simulating several 
sets of Gaussian distributed Xi on the gold set SN positions 
and obtaining A^2 from each realization. For comparison 
with theory we follow GSL08 and compute a bootstrap dis- 
tribution by shuffling the data values Zi, jJ,{zi) and 0-^(2:;) 
over the SNe positions (for further details see GSL08). 

As mentioned in GSL08, the above procedure separates 
the data sets into hot and cold SNe that have large and 
small dispersions with respect to the best fit model. How- 
ever, note that these two sets could still indicate the same 
cosmology, albeit with a different value of x^- To ameliorate 
this deficiency we now introduce a new statistic that does 
not suffer from this artifact. 



3.2 



A^ statistic 



As mentioned above, Xi does not contain information about 
whether the SN is above or below the fit. An obvious gen- 
eralization that does contain information regarding whether 
the SN at a redshift is closer or farther from us can be ob- 
tained by considering a statistic based on XjS. We consider 
two subsets of data defined by two hemispheres labeled by 
the direction vector n, containing A^north and A'aouth SNe, 
where the total number of SNe, A'^ = A'north + A'south, and 
define the quantity 

E 



Axfi 



E 






Xl 



(11) 



Clearly (Axn) = and ((Axn)^) = 1. From the central 
limit theorem (JKendall fc Stuarj U-QTJ ) it follows that for 
N ^ 1, the quantity Ax follows a Gaussian distribution 
with a zero mean and unit variance. As in the previous case 
we maximize this quantity by varying the direction h across 
the sky to obtain the maximum absolute difference 



max{|Axn|} 



(12) 



Unlike the A^2 statistic this statistic is not marginalised 
over the Hubble constant since our results show that 
marginalising over it instead of using its best fit value has 
only a marginal effect on A^2 . Moreover, in the limit N ^ 1, 
and assuming a uniform sky coverage, we expect the two 
hemispheres to contain roughly an equal number of SNe. In 
this case it is clear that A^ would depend only weakly on 
Ho. 

This statistic differs from the previous one in that the 
A;^ statistic has a theoretical limit where the position and 
the shape parameters can be determined analytically. Given 
Nd independent directions on the sky we are essentially de- 
termining the maximum of a sample of size Nd where the 
individual numbers are drawn from a Gaussian distribution 
with a zero mean and uni t variance. In the limit N d 3> 1 the 
parameters are given by (JHaan fc Ferreira Il200a ) 



\/21ogAfd-loglogAfd 
1 



log 47r 



(13) 
(14) 



where we have to additionally assume that the number of 
SNe A'' ^ 1, since the distribution for x becomes Gaussian 
only in this limit. This is convenient since at least for large 
data sets, which will be available in the future, a comparison 
with theory becomes simpler. However, for a smaller num- 
ber of SNe there is a possibility that not all directions are 
independent, in fact, it is quite possible that two directions 
contain exactly same subsets in the two hemisphere. In this 
situation is is clear that the total independent directions is 
a smaller number than Nd and thus theoretical distribution 
would be rightward shifted and also more sharply peaked. 
For this reason we also calculate the bootstrap distribution 
and the theoretical distribution in the same manner as for 
the previous statistic. 



4 RESULTS 

In GSL08 we discussed a specific bias in the bootstrap dis- 
tribution, showing that it is shifted slightly to the left of 
the theoretical distribution due to the fact that theoreti- 
cal distribution is obtained by assuming XiS to be Gaussian 
random variates with a zero mean and unit variance. Theo- 
retical XiS are unbounded, however the bootstrap distribu- 
tion is obtained by shuffling through a specific realization 
of Xi where the x^s are obviously bounded. It is clear that 
on the average this should produce slightly smaller values 
of A-1^,2 in comparison to what one expects from a Gaussian 
distributed XjS. The same problem persists in the Hubble 
constant marginalised version of A^2 . 

However, the analysis in GSL08 had a numerical bug 
that produced a mismatch between the bootstrap and the- 
oretical distributions that is worse than what one obtains 
upon correction. Therefore, for reference we plot the results 
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^2 

Figure 3. The theoretical and the bootstrap probabihty distri- 
butions for GD04 for the A 2 statistic. Theoretical distribution 
is shifted to the left compared to what we find for our simulated 
data in Figure [2] which uses Gaussian deviates suggesting evi- 
dence for non-Gaussianity. 



Table 1. The model parameters for GD04 and GD07 are tabu- 
lated here. 



Model 



Set 



ACDM 
ACDM 



GD04 0.32 
GD07 0.33 



177.1 
158.7 



Table 2. Direction for maximum A GD04 and GD07 are tabu- 
lated here. 



Model 


Set 


A^. 


long 


lat 


ACDM 
ACDM 


GD04 
GD07 


0.83 
0.53 


90. 
347. 


44.7 
27. 



for a new simulation in Figure [5] with a total of 157 SNe. 
This may be compared with Figure 3 of GSL08. The dis- 
crepancy has been reduced after correction. Our results in 
this paper should be interpreted with respect to Figure [2] 
in this paper. Concerns regarding the small number of SNe 
in data sets and its effect on the efflcacy of our method 
can be addressed (as in GSL08) by noting that this figure 
is produced with only 157 SNe and the theoretical and the 
bootstrap distributions look similar. 




Bootstrap 
Theoretical 



0.5 



1.5 



&X^ 



Figure 4. The theoretical and the bootstrap probability distri- 
butions for GD07 for the A 2 statistic. Theoretical distribution 
is compatible with Figure [5] for the simulated data, which implies 
Gaussianity for the residuals. The A 2 for GD07 is within one 
sigma of the mode and thus the data are consistent with isotropy. 



4.1 Results: A^2 Statistic 

This statistic is different from the one used in GSL08 in 
the fact that here we have marginalised over the Hubble pa- 
rameter. However, our results differ significantly from those 
presented in GSL08 due to the fact that we have corrected 
a numerical bug in the calculation used in GSL08 [j In Ta- 
ble [l] we give the best fit values of SIm using the likelihood 
function in Ea. lA7l for both data sets. This is the model that 
we subtract from data to produce the residual XjS. We note 
that GD07 gives slightly higher value of Qm compared to 
GD04. The direction and the value of ( A^2 ) is presented in 
Table [1 

GD04: In FigOwe plot the bootstrap and the theoreti- 
cal distribution expected for GD04 and mark the position of 
GD04. Comparison with Figure 1 in GSL08 shows that for 
GD04 the essential difference is that the theoretical curve 
has shifted to the left while the bootstrap distribution is 
almost identical, suggesting that the effect of marginalisa- 
tion over the Hubble parameter is small. Comparison with 
Figure [T] shows that there is still a weak signature of non- 
Gaussianity since the theoretical distribution instead of be- 
ing to the right of the bootstrap is instead shifted to the 
left. 

GDOT Results are plotted in Fig[4l This should be com- 
pared to Figure 2 of GSL08. As in the case of GD04, we find 



^ In GSL08 the theoretical and the bootstrap distributions were 
handled by different codes, and one of them had a bug thereby 
creating a greater discrepancy between them than should have 
been. 
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that the theoretical distribution has shifted to left after cor- 
recting the bug in the code while the bootstrap distribution 
is unchanged. However, in this case, due to the correction 
of an additional mistake in the code, GD07 has shifted its 
position to the left. Comparison with Figure [2] shows that 
our revised results are compatible with the absence of any 
features of non-Gaussianity in the data. GD07 sits at about 
one sigma from the mode of the distribution, thus the direc- 
tional dependence is weaker than found in GD04. 

4.2 Results: A;^^ statistic 

The essential difference between the previous statistic and 
this one is that this statistic is sensitive to a SN being above 
or below the best fit. As mentioned in § 3.2, in the limit A'^ 3> 
1 and Nd 3> 1, the distribution for A;^ is fully determined. 
In the results below we call this the analytic distribution. 

Since this statistic is not marginalised over the Hubble 
constant, as a first step we fit the data sets to the two pa- 
rameters Ho and f^M in order to obtain the residuals Xi- We 
tabulate the best-fit value of these parameters in Table [S] 
The value of A;^ and the direction in which the maximum 
occurs are tabulated in Table 3] The quoted directions in 
both cases refer to a negative value A;^, thus in both cases 
the direction points in the direction where SNe appear closer 
to us than the best fit. 

GD04: The bootstrap, theoretical and analytic distri- 
butions are plotted in FigO The analytic distribution is very 
different from the other two. However, this can be explained 
by the fact that our assumption of theoretical limit may not 
be satisfied in this case due to a small number of SNe in 
the data. However, the bootstrap distribution agrees quite 
well with the theoretical distribution thus indicating Gaus- 
sianity. Similar to the previous statistic, we find that GD04 
is slightly more than 2a away from the mode of the distri- 
bution. One thing to note is that here A^ is smaller than 
the mode of the bootstrap distribution, indicating a smaller 
anisotropy in the data than is expected from a purely sta- 
tistical point of view. 

GD07: The results are plotted in Fig [6] The main 
difference from GD04 is that the anisotropy is smaller in 
this case. Therefore, GD07 is consistent with Gaussianity 
and isotropy, in agreement with what we find with the A^2 
statistic. 



5 CONCLUSIONS 

We have presented corrected results for the A statistic in- 
troduced in GSL08 and showed that the results for the A^2 
statistic change substantially after correcting for a numerical 
bug in the our code. Our present method extends the previ- 
ous work by invoking an estimator that is marginalised over 
the Hubble parameter. However, our results are not affected 
substantially due to this change since the bootstrap distri- 
bution and the location of A^2 does not change substan- 
tially for GD04 compared to our previous results. However, 
for GD07 we have corrected another mistake in our earlier 
analysis and we find that the bootstrap distribution as well 
as the position of GD07 with respect to it have changed. 
Our main conclusions are 1) the match between the boot- 
strap and the theoretical distribution is much better than 




Figure 5. The theoretical, analytic and the bootstrap probability 
distributions for the A^ statistic for GD04. The analytic distribu- 
tion uses the limiting values of the shape and position parameters 
for the Gumbel distribution. Since the total number of indepen- 
dent directions is most likely smaller than N^ , we overestimate the 
position of peak. The bootstrap distribution indicates anisotropy 
at the level of slightly greater than about two sigma. 




Figure 6. The theoretical, analytic and the bootstrap probabil- 
ity distributions for the A^^. statistic for GD07 as in Figure [5] 
The bootstrap distribution indicates consistency with isotropic 
distribution of SNe in GD07. 
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Table 3. The model parameters for GD04 and GD07 are tabu- 
lated here. 



The inverse transformation is 



Model 



Set 



Hq 



ACDM 
ACDM 



GD04 
GD07 



0.30 
0.33 



64.5 177.3 
63.0 158.8 



Table 4. Direction for maximum A GD04 and GD07 are tabu- 
lated here. 



Model 



Set 



long lat 



ACDM 
ACDM 



GD04 
GD07 



-1.52 

-1.77 



107 
97 



26 

70 



was presented in GSL08 2) GD04 shows some evidence for 
non-Gaussianity, however, GD07 is statistically consistent 
with a Gaussian distribution of residuals 3) and although 
GD04 has a weak direction dependence, GD07 is consistent 
with an isotropic distribution of SNe. 

Our new statistic A^ similarly shows a weak direction 
dependence in GD04 but no significant anisotropy in GD07. 
One surprising feature we find is that for both GD04 and 
GD07, A-,^2 turns to be larger than the mode of the distribu- 
tion, while Aj(. is smaller than the mode of the distribution. 
However, since our results are consistent with isotropy, we 
do not investigate the implications of this puzzling feature. 
T he results of this statistic ca n be compared to the results 
of lCooke fc Lvnden-Belll ([20101) . The directions of maximum 
anisotropy are different. However, since the methods and 
data are different, and the results suggest only a weak di- 
rection dependence, this mismatch is not unexpected. 



APPENDIX A: MARGINALISATION 

The normalized likelihood function for the distance modulus 
is given by 



L(x) = 



(27r)^/2x/detS 



exp 



2^' ^ 



(Al) 



Here Xi = fii — fJ,thcoiy{zi; Ho, Om, f^x), where the subscript 
j = 0, . . . , A'' and E^j = Sij a'f is the covariance matrix. The 
superscript 'T' denotes the matrix operation of taking the 
transpose. Each modulus has a linear constant term con- 
taining a combination of parameters Hq and M, therefore, 
we wish to consider the difference of two moduli to explic- 
itly remove it. For concreteness we choose the O"' SN as 
the anchor point, subtract its modulus from all others, and 
marginalise over it. We first define these new parameters y 
as follows. 



y = A-^ 

where the transformation matrix 
1 



(A2) 




for i — j 

for i > 0, j = 

otherwise 



x = Ay. 



(A3) 



Noting that det A = 1, in terms of the new variables, the 
likelihood function is given by 



L{y) = 



(27r)'v/2\/detS 



exp 



1 To-l 

-^y s y 



(A4) 



where S"^ = A'^S'^A. 

We now integrate over yo to obtain the marginalised 
likelihood function. 



L{yi,y2,...,yN) 



1 / 1 T^-l 

, exp V C y 

(27r)^/2%/ditC V 2-^ 



where, the final covariance matrix C is given by 

-1 "^Oi '^Oj 



CT^^SZ 



^00 



(A5) 



(A6) 



where the indices i,j run from 1 to A'^. To formulate 
our estimator for H^ marginalised statistic we require the 
marginalised Xm defined as follows 



2 T^-1 

Xm = y C y 



(A7) 



We note that the expectation value of (xl/) ~ N as shown 
below 



ixli) = (y^C- V) = tr(C-^C) = N 



(A8) 



In this formulation the marginalisation is carried out using 
one of the SN data points. However, it is easy to see that 
the method is more general and we can marginalise over any 
Gaussian prior on the Hubble parameter. To achieve this 
we note that at low redshifts the luminosity distance takes 
the approximate form cLl — cz/Hq. Thus, it is possible to 
produce a theoretical luminosity distance at low redshifts 
with appropriate error bars if we are given Gaussian priors 
on the Hubble parameter and the absolute magnitude for 
type la SNe. If we use this theoretically produced magnitude 
and replace the O"' SN with it then we can marginalise over 
any arbitrary priors on the Hubble parameter. 
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